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We have studied the correlation effects in Cs and Fr arising from the interplay of the residual 
Coulomb interaction to all orders and the neutral weak interaction which gives rise to the parity 
violating electric dipole transition to first order, within the framework of the relativistic coupled- 
cluster theory which circumvents the constrain of explicitly summing over the intermediate states. 
We observe that, the contributions arising from the perturbed doubly excited states are quite signif- 
icant and hence, any calculation should not be considered accurate unless it includes the perturbed 
double excitations comprehensively. In this article, we have reported a comparative study of various 
results related to the parity violation in Cs and Fr. 
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I. INTRODUCTION 

One of the most challenging contemporary problems in physics is the search for possible new physics beyond the 
well known standard model (SM) of elementary particles [l|, H, H, 0|- Apart from using gigantic accelerators at 
high-energy scales, it is also possible to use high precision, albeit, low-energy table top atomic experiments, such 
as, the measurement of atomic parity violation, in combination with accurate relativistic many-body calculations 
of the atomic parity nonconserving transition amplitude, Ei PNC to achieve this goal 0, @- Some of the prominent 
signatures of physics beyond the SM which can be inferred from these atomic experiments are: a tight limit on the 
mass of extra Z-bosons^precise value of the Weinberg angle, limit to the radiative corrections for the electron-nuclear 
weak interactions etc [J]. As some of these SM results are known to high precision, they demand similar sub-one 
percent accuracies in both the measurements and the atomic calculations. As the interaction Hamiltonian for the 
atomic parity violation (APV) due to the nuclear spin-independent (NSI) electron-nucleus interactions is proportional 
to Z 3 where Z being the atomic number Q, heavy atomic systems are chosen for the study of APV effects. A series of 
APV experiments on a number of atomic systems, including those of Cs @ and Tl 0] have been carried out. However, 
the high accuracy of ~ 0.35% has been achieved only for Cs. Furthermore, a number of ab initio calculations of APV 
amplitudes in Cs have been carried out using a variety of many-body approaches. Some results based on the relativistic 
coupled-cluster (RCC) theory are also available for Cs, however, their accuracies are somewhat uncertain, since most 
of them have used the sum-over-states approach which considers contributions from the core orbitals approximately 
and accounts for only a selected number of excited states whose contributions are dominant. In addition, the doubly 
excited intermediate atomic states and the normalization of the RCC wave functions are treated only approximately. 
We have developed a technique in the frame work of RCC theory that circumvents these drawbacks and it has been 
employed earlier in the calculation of the EIpnc amplitudes in Ba + Q and Ra + in which we have demonstrated 
that the accuracies of < 1% and < 3%, respectively, were possible. 

In this work, we employ the new RCC approach, mentioned above, to study various correlation effects in the parity 
violating amplitudes in Cs and Fr. We report a comparative study of their results along with those reported previously. 

II. THEORY OF APV 

The dominant interaction in an atom is the electromagnetic interaction which, as well known, conserves parity. 
However, there is a non-zero probability of the interaction between the electrons and the nucleus of an atom due to 
the weak force with the exchange of a Zq boson, as shown in Fig. [TJ which violates parity. The interaction Hamiltonian 
between the electrons and the nucleus due to the weak interaction can have two components: one, vector-axial- vector 
and the other, axial-vector-vector currents. The latter depends on the nuclear spin and most of its contribution 
cancels out except from the odd nucleon and hence, it is relatively smaller in magnitude [H than the former; the NSI 
component. In this work, we shall consider the APV effect due to the NSI component alone. 
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FIG. 1: Diagrammatic representation of the electron-nucleus interactions due to the electromagnetic (with the exchange of 
photon 7 denoted by solid curved line) and the weak (with the exchange of heavy boson Zo denoted by straight dashed line) 
interactions. 



The APV interaction Hamiltonian due to the NSI component is given by, 

H^v = ^=Q w (Nh!p%(r), (2.1) 

where Gf is the Fermi constant, p e N (r) is the electron density over the nucleus and 7§(= *7o7i727f)i is the product 
of the four Dirac matrices that involve electron spin, and Q W (N) is the nuclear weak charge, which is equal to 
2(Zcp + Ned) where c p and c„ denote the electron-proton and the electron- neutron coupling constants for the atomic 
(Z) and neutron (TV) numbers, respectively. The values of these coupling constants predicted by the SM, in the lowest 
order of electroweak interaction (at the tree level), are given by 

c p = \{l - 4sin 2 6 w ) a 0.04, c„ = ~, (2.2) 

where 6\y is the Weinberg angle and its measured value is ~ sin 2 8\y w 0.23 [Toj . Substituting these values in Q W (N), 
we get Q W (N) = — N + Z(l — sin 2 9 w ) which is proportional to N. Hence, the perturbation due to H^p V is generally 
expressed in the scale of Q W (N)/N. Since H^p V does not commute with the parity operator, its inclusion with the 
atomic Hamiltonian of the electromagnetic interaction, which commutes with the parity operator, mixes the opposite 
parity states of same angular momentum. The strength of this interaction is sufficiently weak, which justifies its 
consideration as a first-order perturbation. 

The Dirac-Coulomb (DC) Hamiltonian used, here, in the calculation of the atomic wave functions of definite parity 
is given by 



H DC = [cat ■ Pi + l)c 2 + V nuc (n)] 



i>j 



(2.3) 



where c is the velocity of light, a and (3 are the Dirac matrices (note that j e — [3a) and V nuc (r) is the nuclear 
potential. 

The atomic wave functions (l^l ^)) corresponding to H DC can be considered as the unperturbed wave functions. 
The total wave function of a system including the first order correction due to the interaction Hamiltonian H^py is 
given by 

\* n ) = \*M)+G F \¥V), (2.4) 

where is the first order perturbed wave function of its unperturbed valence state | V I'1 ' > ) and Gf is used as a 

coupling constant. 

Electric dipole (El) transitions between the states of same parity are forbidden due to the electromagnetic selection 
rules. However, an El transition between the states of mixed parity, mixed due to H^py interaction, is possible and 
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the corresponding transition amplitude can be expressed as 

E1 = (* f \Dm 

PNC V(*/I*/X*il*i>' 

where D = e r is the El operator, the subscripts i and / denote initial and final valence orbitals, respectively. 

Expanding the total wave function as given in Eq. (|2.4p and retaining the terms only up to first order in Gp, we 
get 

(*f)\D\^) + {^\D\^) 
rji-PNC = Gp- 



■*(°)|$(p))(*(°)|$( )) 



(^ 0) |J)|^ 0) )(^ 0) |g^|^ 0) ) 
Ei — Ei 



+Gf 2^ eT^Ej ' (2 - 6) 

J=£f J 

where the subscripts I and J represent the intermediate unperturbed states. We have used, here, the explicit form 
for the first order wave function given by 

An important question we address in this paper is: How significant are the contributions from those states which 
were considered approximately in the sum-over-states approach and how they vary with the size of the systems? We 
would address this by carrying out a comparative study of EIpnc results in two systems, namely Cs and Fr, of 
different atomic sizes. In order for the contributions of these higher excited states to be included, it is necessary to 
solve the first order perturbation equation directly. In other words, it is necessary to solve the equation following 
equation 

{H oc _ £(0)^(1)) = x_ {E p> _ H$&)\*W), (2.8) 
where E^(= \H/[pv I*" '')) i s the first order correction to E^ which, however, vanishes in the present case. 



III. APPLICATION OF RCC THEORY TO APV 



The RCC method, which is equivalent to all order perturbation theory, has been used in the recent past and accurate 
results have been reported for many single valence systems [1, S [0, H2] • I n the RCC framework, the wave function 
of a single valence atom can be expressed as 

= =e T<0) {l + ^)}|$„), (3.1) 

where |$ n ) is the reference state constructed from the Dirac-Fock (DF) wave function |4>o) of the closed-shell config- 
uration by appending the valence electron n, that is, |$„) = a/ n |3?o) where represents a creation operator which 

creates the valence electron n. Here and Sn^ are the RCC excitation operators which excite electrons from 
|<E>o) and \&n), respectively, due to the residual Coulomb interactions. The corresponding excitation amplitudes are 
obtained by solving the following equations 

($ L \{H™^ W }\4> ) = (3.2) 
(^\{H^ w }S n \^ n ) = -(^\{H^°^\^ n ) + (^SMAE^, (3.3) 

with the superscript L(— 1, 2) representing the singly and doubly excited states from the corresponding reference 
states and the wide-hat symbol over H® c e Ti0) represent the linked terms of normal order atomic Hamiltonian H^ c 
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and RCC operator T^°\ In the CCSD (CC with single and double excitations) approximation, the corresponding 
RCC operators are defined by 

T (o) = r (0) + T (o) (3 4) 

and 

for the closed-shell and single valence open-shell systems, respectively. The quantity AEn^ in the above expression 
is the electron affinity energy (or negative of the ionization potential (IP)) for the valence electron which is evaluated 
by 

A£f = {* n \{Hjfi? m }{l + S®}\* n ). (3.6) 

In addition to having considered full singles and doubles in the CCSD equations given in Eq. (|3 . 3[) . we have also 
included the contributions from the important triple excitations perturbativcly (known in the literature as CCSD(T) 
method) by defining 

s (V, P ert = H D^O) +H -j^» j (37) 

where the superscript pert denotes the perturbation and their contributions to AEn^ are evaluating as 

AE (0),tri P = T Wt~^W, P ert (3.3) 

After solving for the amplitudes of T™, we solve Eqs. (|3.3|) and ()3.6() simultaneously and obtain the amplitudes of 
Sn* 1 operators. 

Now, the total atomic wave function in the presence of H^py 1S expressed, in the RCC ansatz, as 



*n> = e TW + G ^ w {l + S^+G F S^ n ), (3.9) 



where and S-h, are the first order perturbed amplitudes corresponding to the unperturbed RCC operators 
and Sn , respectively. On expanding the above equation keeping the terms only up to first order in Gp yields 

= e T(0) {i + s^ + G F T^{i + s^) + G F s^}\^ n ). (3.10) 

Comparing the above equation with Eq. (12. 4p . we get 

|*W> = e T<0) {T«(l + ^)) + ^}|$„). (3.11) 

In order to calculate \^n) as a solution of Eq. (|2.8p in the RCC theory, we solve the excitation operator amplitudes 
of and S { n 1] using the following equations 

(<f> L \{H^ m T^}\<S> ) = -($Vll^ <0, |<i>o) (3.12) 

and 

(^\{H^ m }S^ n ) = -(^\{H^ m T^(l + S^) 

+H^e~T m (l + S< >)}|$ n ) + (^\SP\^ n )AE n , (3.13) 

after solving Eq. (|3.2p and Eq. (|3.3p . respectively. In the above expression, notation H^pye T(0) is used for the 
connecting terms between Hj[py and T^°\ To keep the level of approximation uniform through out, both T"W and 
Sip are truncated at single and double excitations by defining 

T (i) _ T (i) +T (i) ^ 314 ^ 



5 



and 

S« = (3.15) 

where and correspond to the perturbed single excitations and T^p and S^J correspond to the perturbed 
double excitations, from closed- and open-shells, respectively. Since both the perturbed single and double excitation 
amplitudes are solved simultaneously, certain correlation effects due to the perturbed double excitations also reflect 
indirectly in the contributions of the perturbed single excitations. 

After obtaining the unperturbed and the perturbed RCC operator amplitudes in both the closed-shell and one- 
valence open-shell atoms, we proceed to calculate the EIpnc amplitude as 



EIpnc = Gi 



= q (gnjjl + S^}D{T^\l + Sii") + S£}|<fr n ) 
<$„|{l + 5i 0)t }7V {l + 5i 0) }|<I> n ) 

+&F , nH _ jm , (3.16) 

<$„|{l + ^ 0)t }7V {l + ^ 0) }|$ n > 

where we define D = (e T< ° )+ Z?e T< ' ) and TVo = e T<0)t e T< °' . The non-truncative series for D and No are expanded using 
the Wick's generalized theorem and are truncated when the terms are below fifth order of the Coulomb interaction. The 

core- valence and valence correlation contributions are obtained from Dt[ 1] andD(T 1 (1) {l + S*i 0) } + 5^ ) ), respectively, 
along with their conjugate terms. 

Corrections due to the normalization of the wave functions are accounted by evaluating 



Norm 



{^\D\^) + {^\D\^) 



N„ 



1 + N n 



(3.17) 



where N n = ($„|{1 + S^}N {1 + Si 0) }\$ n ). 

Although, the CCSD(T) method described here, accounts for the contributions from the important unperturbed 
triple excitations it fails to include the direct triple excitation contributions to the EIpnc calculations. To account for, 
at least, the lowest order direct triple excitation contributions (minimum up to fourth order in Coulomb interaction), 
we construct them with the open-shell RCC operators perturbatively as follows 

, , zjDCt(O) i TjDCqiO) 
b nab - ~ ~ ~ 1 Z T> (S.lB) 



+ ( a + £b — £p — £q 



and 



, N ttDCt (1 ) 4- TjDCqi 1 ) 

Q pqr,(l) _ -"jy ± 2 1-n N o 2n 

nab ~ c J- C J. <=, - c - c - c ' ^ ' 

"r to i t-b fcp fcr 

where is the single particle energy of an orbital i. These operators are finally considered as parts of Sn^ and in 



Eq. piip . This can be called as lo-CCSDvT approximation, which also accounts important lower order valence triple 
excitation effects in the final property calculations. 

IV. RESULTS AND DISCUSSIONS 
A. Orbitals generation 

We have used Gaussian type functions 

F GTO {r t ) =r n «e- Qir . ? (4.1) 
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to construct the DF orbitals where ccj is an arbitrary parameter which has to be chosen and represents a radial 
grid given by 



n = r 



,h(i-l) _ l 



(4.2) 



where the step size h is taken to be 0.03, the radial grid is increased up to i = 750, tq is 2 x 10 6 in atomic units and 
n K is the radial quantum number of the orbitals. Here aiS are chosen to satisfy the even tempering condition 

a, = ao/3 1 " 1 (4.3) 

and we have chosen different ao and (3 values for different symmetries (I) (known as even tempered basis) and they 
are given in Table [U 

TABLE I: Used qq and f3 parameters for different symmetries (I) to construct DF orbitals using GTOs in Cs and Fr. 



ao 



0.00190 0.001825 0.00183 0.00185 0.00187 
2.91 2.94 2.94 3.05 3.09 



The finite size of the nucleus in these systems is accounted by assuming a two-parameter Fermi-nuclear-charge 
distribution for evaluating the electron density over the nucleus as given by 

p( n ) = f° , (4.4) 

' y tJ i _|_ e (n-c)/o' ^ > 

where po is the density for the point nuclei, c and a are the half-charge radius and skin thickness of the nucleus. These 
parameters are chosen as 

a = 2.3/4(Zn3) (4.5) 

and 



C =\j\ r rrns-\°?K 2 , ( 4 - 6 ) 

where r rms is the root mean square radius of the corresponding nuclei which is determined by 

r rms = 0.836A 1/3 + 0.570 (4.7) 

in fm for the atomic mass A. 

B. Results 

In Table [Til we give the Elpjyc amplitude results obtained from various calculations. From this work, we present 
the results using the DF, CCSD, CCSD(T) and lo-CCSDvT methods for both Cs and Fr. In the same table, we 
also compare our results with previously reported results using various many-body methods. Our Cs result matches 
reasonably well with the other calculations, but our Fr result differs significantly. We give below individual contribu- 
tions from various RCC terms and express them in terms of correlation diagrams and level of excitations in order to 
facilitate the readers to understand the role of various correlation effects and intermediate states. Briefly, we discuss 
here the methods used in the other calculations. The difference between our DF and lo-CCSDvT results gives an 
idea about the amount of total correlation effects through the present method in these calculations. Considering the 
same DF wave functions, we have also employed CCSD and CCSD(T) methods where we find that the CCSD results 
are larger in magnitude than the CCSD(T) results. However, the EIpnc result increases in magnitude for Cs in 
the lo-CCSDvT approximation, but it decreases for Fr; indicating that the triple excitation effects are stronger in Fr. 
From a comparison between these results, we observe that the dominant triple excitations effects arise through the 
CCSD(T) method. 
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TABLE II: Comparison of EIpnc results of Cs and Fr due to electron-electron Coulomb interactions from various methods 
in xW~ n ieao(-Qw/N). 



Cs 223 Fr Method Reference 



'1/2 — ► 7s 5*1/2 


7s 2 Si /2 8s 2 5i /2 






-0.7236 


-13.4825 


DF 


This work 


-0.8941 


-14.5134 


CCSD 


This work 


-0.8889 


-14.4219 


CCSD(T) 


This work 


-0.8892 


-14.4106 


lo-CCSDvT 


This work 


-0.904(9) 




LCCSD(T) + SL + SS 


[13] 


-0.9078 




BO + GFCP 


[14] 


-0.904 


-15.72 


CI 


[15] 


-0.902(4) 




CCSD(T) 


[16] 


-0.9001 


-15.229 


BO + GFCP 


[17] 


-0.8998(25) 




CCSDvT + SS 


[18] 




-15.9 


BO + GFCP 


[19] 




-15.41(17) 


LCCSD + RPA + ExpEng + SS 





NOTE: Contributions from Breit interaction, QED correction and nuclear effects are not considered here. 

Abbreviations SS : sum-over-states 

SL : scaling 

BO : Brueckner orbitals 

GFCP : Green function technique for all order correlation potential 

CI : configuration interaction method 

ExpEng : experimental energy 

RPA : random phase approximation 



C. Discussions 



Here we discuss briefly different reported calculations of the above results at the DC approximation. About two 
decades ago, Blundell et al [l3| had employed the linearized CCSD(T) method (LCCSD (T) method) to evaluate 
the unperturbed wave functions in Cs and then they had used a sum-over-states approach to evaluate the Elp^c 
amplitude for the 6s 2 S 1 / 2 — > 7s 2 S , 1 / 2 transition in Cs. However, they had scaled their wave functions to fit the 
calculated energies of different states with the experimental results which reproduced many atomic properties quite 
accurately, but that does not show that the method they had used is capable of producing accurate ab initio results. 
Shabaev et al [l5| have also obtained results using a CI method with a local form of the DF wave functions. Their Cs 
result matches with the result of Blundell ct al. In our previous work on Cs |l6[, we had also calculated this quantity 
by considering the same «o and (3 parameters for orbitals of all symmetries (known as universal basis). With the 
new parameters, the Gaussian basis orbitals produce better wave functions in the nuclear region which are verified 
by studying the hyperfine interactions that will be reported elsewhere. The convergence of the RCC amplitudes are 
better than 10 -8 in the present case than 10~ 6 in [l6[ and it gives a slightly different result. Dzuba et el have carried 
out a few calculations of these quantities using Brueckner orbitals using a Green function technique (Feynman diagram 
approach) that takes into account various classes of correlation effects to all orders and avoids the sum-over-states 
approach 0, [l?], Their results also differ from each other and in some cases with others as can be seen in 

Table fl] The most recent calculation on Cs is reported by Porsev et al [l8[ using the RCC method that includes 
all single and double excitations with all valence triple excitations (CCSDvT method). However, they have finally 
used a sum-over-states approach to calculate the Elp^c amplitude of the 6s 2 S*i/2 — > 7s 2 Si/ 2 transition. They give 
contributions from 6pi/ 2 to 9pi/ 2 singly excited states as "Main", which contributes —0.8823(18), and the remaining 
contributions as "Tail", which is obtained using other many-body methods as —0.0175(18), (results are always given 
in xl0~ n ieao(— Qw/N) here onwards) at the DC approximation. In contrast to their approach, we have used the 
lo-CCSDvT method, but have included contributions from the core and doubly excited states in a manner similar to 
that of the singly excited states. 
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TABLE III: Contributions to the EIpnc calculations in Cs and Fr using CCSD(T) method (in xlO" n iea (-Qw/-/V))- 



(A) Cs 
Initial pert 



terms 

2 c (0) ■ 



Final pert. 



terms 



DF result 

T-irriVSi 



-1.0460 



0.3224 



T^D 0.0374 DT (1) -0.0355 

DiS^J -1.7467 S^Di 0.1841 

^i^P -0.0046 SffiDi 0.0226 

S^tOiS^ 0.4061 S^D-lS^ 0.2141 

S'^ )t DiS'^ ) -0.0055 S^SiS^ 0.0301 

S^DiS^ -0.0026 S^UiS^ -0.0028 

S^ f DiS^ } 0.0011 S^SiS* -0.0073 

Ot/iers " -0.0045 " 0.0054 

Norm 0.0301 -0.0105 



(B) Fr 






Initial pert, terms 


Final pert. 


terms 




(7s>S[% - 




DF result 






UttAPV 


-19.4661 H^ V D 


5.9836 



T^D 0.7147 DT {1) -0.6876 

DiS[V -31.7107 S^Ui 4.0538 

TliS^ -0.0314 5^ 3t Si 0.5930 

SWDiS^ 7.8001 S^DiS^ 4.1072 

S&DiSiV -0.0423 S^UiS^ 0.6580 

S^DiS^ -0.0615 S^SiS™ -0.0911 

S^DtS^ 0.0033 S^SiS^ 3 -0.1501 

Ot/iers " -0.0152 " 0.0769 

Norm 0.6043 -0.2433 



D. Cs vs. Fr 



In Table IIII| we present the individual contributions from different RCC terms in our CCSD(T) method to the 
EIpnc calculations for the 6s 2 Si/ 2 — ¥ 7s 2 5i/ 2 an d 7s 2 S 1 / 2 — > 8s 2 Si/ 2 transitions in Cs and Fr, respectively. It 
is evident that the trends of the contributions from the different RCC terms are similar for both Cs and Fr. The 
important contributions come from three terms: DiS[]\ D%S^' and DT[ x ^ and their corresponding conjugate terms 
where D\ is the effective one-body terms of D and bare operator D is its lowest order term. Other terms correspond 
to higher order RCC terms, but they are not small. Contributions given as Others represent the RCC terms that 
come from the effective two-body terms of D after contracting with the open-shell RCC terms. We give diagrammatic 
representations of the above three RCC terms in Fig [5] (without their conjugate terms) and their lowest order terms. 
From this figure, it can be noticed that DT[^ and DxS[l contain the lowest order DF contributions from the core 

(hole) and virtual orbitals, respectively; hence DiS$ always gives the largest contribution. Again, the perturbed 
states arising through the ground state contribute predominantly while contributions from the perturbed excited 
states are comparatively smaller than those corresponding to the ground state but with the opposite signs. The 
final results are the outcome of these cancellations. The most core correlation effects are coming from DT^ and 

its conjugate terms while there are small contributions that come from DT^ and its conjugate terms which are 
included in "Others". As seen these correlation effects mostly cancel out with their corresponding conjugate terms in 
both the systems. The contributions from singly excited states can be estimated by summing the total contributions 
from DiS[]\ s[°^ DiS[]\ sffi DiS-j^ and their conjugate RCC terms. The remaining contributions are come from 
the doubly excited states apart from the normalization corrections to the RCC wave functions. After accounting the 
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TABLE IV: Important contributions to the EIpmc calculations in Cs from various intermediate states (I and J notations are 
used as per Eq. (|2.6p in x 10 _11 ieao( — Qw/N)). 





-> is o 1/2 


I 


6s*S\% - 

Results 


-> /s o 1/2 


J 


Results 





0005 




Ari-i i~ — > fie 


—0 0005 


5Pl/2 — * 7s 


0.0368 




5pi/2 — » 6s 


-0.0349 




S^D, 




6s -> 6pi /2 


-1.9195 




7s -> 6pi/ 2 


1.7648 


6s — > 7pi/ 2 


0.1532 




7s — > 7pi/ 2 


-1.5022 


6s -> 8pi /2 


0.0347 




7s -> 8pi /2 


-0.0902 


6s -> 9pi /2 


-0.0084 




7s -> 10p 1/2 


0.0093 


6s -> 10pi/2 


-0.0051 
















6s -> 6pi /2 


0.1716 




7s -> 6pi /2 


0.2556 


6s — > 7pi/2 


0.1302 




7s -> 7pi/ 2 


0.1213 


6s -> 8pi /2 


0.0258 




7s -> 8pi /2 


-0.0611 


6s -> 9pi /2 


0.0735 




7s -> 9pi /2 


-0.0946 


6s -» lOp^ 


0.0053 




7s -> 10p 1/2 


-0.0055 



6s5px/ 2 — > 7s6s 


-0.0019 


7s5p 3/ 


2 — > 6s9s 


0.0002 


6s5pi/ 2 — > 7s8s 


-0.0001 


7s5p 3/; 


j -> 6sl0s 


0.0009 


6s5p 3 / 2 — > 7s8s 


0.0001 


7s5pi /2 


-> 6s6d 3/2 


0.0001 


6s5pi/ 2 — > 7s9s 


-0.0009 


7s5p 3/2 


— > 6s6d 3 / 2 


-0.0003 


6s5p 3 / 2 — > 7s9s 


0.0009 


7s5p 3/2 


-> 6s7d 3/2 


-0.0001 


6s5pi/ 2 — ► 7sl0s 


-0.0014 


7s5pi /2 


-> 6s8d 3/2 


0.0006 


6s5p 3 / 2 — * 7sl0s 


0.0001 


7s5p 3/2 


-> 6s8d 3/2 


-0.0003 


6s5s — ► 7s9pi/ 2 


0.0001 


7s5pi /2 


-> 6s9d 3/2 


0.0019 


6s5s —* 7sl0pi/ 2 


0.0002 


7s5p 3/2 


-> 6s9d 3/2 


0.0003 


6s5pi /2 -> 7s6d 3/2 


0.0006 


7s5pi /2 


-> 6sl0d 3 / 2 


0.0005 


6s5p 3/2 -> 7s6d 3/2 


0.0005 


7s5p 3/2 


-> 6sl0d 3 / 2 


0.0001 


6s5pi /2 -> 7s7d 3/2 


0.0003 


7s5p 3/2 


-> 6s5d 5/2 


0.0044 


6s5p 3/2 -> 7s7d 3/2 


0.0002 


7s5p 3/2 


— > 6s6d 5 / 2 


0.0018 


6s5pi /2 -> 7s8d 3/2 


0.0013 


7s5p 3/2 


-> 6s7d 5 / 2 


0.0008 


6s5p 3/2 -> 7s8d 3/2 


0.0011 


7s5p 3/2 


— ► 6s8d 5/2 


0.0041 


6s5pi /2 -> 7s9d 3/2 


0.0006 


7s5p 3/2 


—> 6s9d 5/2 


0.0048 


6s5p 3/2 -> 7s9d 3/2 


0.0005 


7s5p 3/2 


-> 6sl0d 5 / 2 


0.0003 


6s5pi /2 — > 7sl0d 3 / 2 


-0.0004 


7s4d 3/2 


-> 6s8/ 5/2 


0.0003 


6s5p 3 / 2 — > 7sl0d 3/2 


-0.0002 


7s4d 3/2 


-> 6s9/ 5/2 


0.0005 


6s5p 3/2 -> 7s8d 5/2 


-0.0003 


7s4d 5/2 


-> 6s8/ 7/2 


0.0004 


6s5p 3/2 -> 7s9d 5/2 


-0.0013 


7s4d 5/2 


-> 6s8/ 7/2 


0.0004 


6s5p 3/2 -> 7sl0d 5 / 2 


-0.0012 


7s4d 5/2 


-» 6s9/ 7/2 


0.0009 


6s4d 3/2 -> 7s8/ 5/2 


-0.0001 








6s4d 3/2 -> 7s9/ 5/2 


-0.0004 








6s4d 5/2 -> 7s8/ 7/2 


-0.0005 









corresponding corrections from the normalization of the wave functions to various RCC operators, we get contributions 
as 0.0018, —0.8980 and 0.0073 from the core correlation, singly excited states and doubly excited states, respectively, in 
Cs. Similarly, we get the contributions as 0.0271, —14.7551 and 0.3158 from the core correlation, singly excited states 
and doubly excited states, respectively, in Fr. Clearly, the core correlation and doubly excited states contributions in 
Fr are large. In Cs and Fr, the doubly excited states contributions are around 1% and 2%, respectively, with opposite 
signs. Hence, these contributions should also be accounted accurately for high precision results. 

To understand more clearly about the role of different excited states that play important roles in the EIpmc result 
of Cs through our RCC approach which have also been studied by other methods, we consider some of the important 
RCC terms and analyze their contributions through various excitations level. These results are given explicitly in 
Table ITVl As noticed, only 4pi/ 2 and 5p\/2 core orbitals are important to account the core correlation effects; this 
is because the single particle energy levels between these orbitals with the valence electrons are small. The singly 
excited states from 6pi/2 to 9px/2 contribute predominantly, but there are also significant contributions from some of 
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(c) 

FIG. 2: Break-down of the perturbed RCC diagrams into the lower-order many-body theory (MBPT) diagrams. Here i and 
/ represent the initial and final orbitals, P is the PNC operator and V is the two-body part of the Coulomb interaction 
Hamiltonian. Lines going up and down with single arrows represent virtual and hole (core) orbitals and lines with double 
arrows represent attached valence orbitals to the closed-shell DF wave function. 




(g) (h) (i) 



FIG. 3: Some of the important contributing diagrams from the triple excitations which contribute significantly to the lo- 
CCSDvT method along with their conjugate terms. 

the continuum states. The continuum orbital, 10pi/2, whose density in the nuclear region is large. There are various 
doubly excited states that play important roles, mainly from the 6skp to 7s\d and Iskp to 6sld excitations for any 
arbitrary principal quantum numbers k and 1, in calculating the above quantities but mostly they cancel each other. 
Therefore, a suitable choice of basis functions and an accurate many-body method are necessary for accounting these 
contributions precisely. 

To compare the role of the correlation effects in calculating E\ PNC amplitudes in Cs and Fr, we compare the DF 
results and then their final results. It is about 20% in Cs, but it is only about 7% in Fr. The small correlation 
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effects in Fr is because of strong cancellations between the contributions from the initial and final perturbed states. 
Again, it is interesting to note that the CCSD results were larger than the CCSD(T) results suggesting that there 
were cancellations from the triple excitations. Using the lo-CCSDvT method, these results increase in Cs slightly 
but cancel out in Fr. The triple excitation effects in Fr are strong enough to be taken seriously into account for high 
precision results. In Fig. [3l we give some of the important RCC diagrams (without their conjugate terms) that give 
large contributions through the lo-CCSDvT method. In the future, we would like to study other relevant properties 
in order to determine the accuracy of our calculated EIpnc amplitudes in these systems. 

V. CONCLUSION 

We have applied the relativistic coupled-cluster method to calculate the atomic parity violating effects in Cs and 
Fr which involve the interplay of the long range electrostatic interactions and the short range weak interactions. The 
present approach considers the Coulomb interaction up to all orders and the weak interaction between the nucleus and 
the electrons up to first order. From the detailed analysis of different contributions, it is clear that the doubly excited 
states play an important role in obtaining precise results. These contributions are quite large in magnitude, however, 
with opposite signs, than those of the singly excited states in Fr, resulting in a large cancellation in the final results. 
They are about 1% in Cs, for which the experimental accuracies have been claimed to be around 0.35%. This suggests 
that, it requires a method like our present approach to consider them accurately for high precision calculations. 
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